Rheology of Active-Particle Suspensions 
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We study the interplay of activity, order and flow through a set of coarse-grained equations gov- 
erning the hydrodynamic velocity, concentration and stress fields in a suspension of active, energy- 
dissipating particles. We make several predictions for the rheology of such systems, which can be 
tested on bacterial suspensions, cell extracts with motors and filaments, or artificial machines in a 
fluid. The phenomena of cytoplasmic streaming, elastotaxis and active mechanosensing find natural 
explanations within our model. 

PACS numbers: 87.16. Ac, 87.15.Ya, 87.10. +c 



An active particle [lj, |2| absorbs energy from its sur- 
roundings or from an internal fuel tank and dissi- 
pates it in the process of carrying out internal move- 
ments usually resulting in translatory or rotary mo- 
tion. This broad definition includes macroscopic ma- 
chines and organisms, living cells, and their compo- 
nents such as actin- myosin and ion pumps In 
this paper, we consider the interplay of activity, or- 
der and flow via coarse-grained equations governing 
the hydrodynamic velocity, concentration and stress 
fields in a suspension containing active particles of 
linear size £, at concentration <p, each particle ex- 
erting a typical force / on the ambient fluid, with 
the activity of an individual particle correlated over 
a time tq (say the 'run' time of a bacterium), and 
collective fluctuations in the activity correlated over 
length scales £ and timescales r . Rather than fo- 
cussing on ordered phases , instabilities 0, , or 
patterns (asters, vortices, spirals) formed in such as- 
semblies |6j which our equations are of course ca- 
pable of predicting, we apply them in the isotropic 
phase, with a view to understanding how a system 
such as a biological cell, composed of active elements, 
responds to deformation or mechanical stress. In ad- 
dition to throwing light on full-cell rheometry 0, > 
our equations form the framework for an analysis 
of any experiment probing the mechanical conse- 
quences of biological activity. 

Our simple model makes rather interesting predic- 
tions: An orientationally ordered state of active par- 
ticles has a nonzero, macroscopic, anisotropic stress 
in contrast with thermal equilibrium nematics. Ac- 
tivity contributes an amount 5rj ~ ftcor to the vis- 
cosity, with a sign determined by the type of active 



particle, and always enhances the apparent (noise) 
temperature. The latter greatly enhances the am- 
plitude of the t~ d l 2 long-time tails in the veloc- 
ity autocorrelation. On approaching an orientation- 
ally ordered state, active suspensions with br\ > 
behave like passive systems near translational freez- 
ing, showing strong shear thickening and Maxwell- 
like viscoelasticity. Nonlinear fluctuation corrections 
give a dynamic modulus G*(uj) ~ \fiuj for lu t _1 
observable over a large dynamic range, since r is 
large. Cytoplasmic streaming m which mate- 
rial flows from the depolymerising trailing edge to 
the polymerising leading edge of a crawling amoe- 
boid cell, finds a natural explanation in our model, 
as do elastotaxis and active mechanosensing ^| > 
where cells orient their motion along preferred axes 
of the ambient medium. 

These results follow from equations of motion based 
simply on the conservation law dtg = — V • a for the 
total (particles + fluid) momentum density g(x, t) 
for an incompressible suspension. The stress ten- 
sor er must in turn be determined by constitutive 
relations which emerge from an additional equation 
of motion for an active order parameter field. We 
ignore, for simplicity, the dynamics of the active- 
particle concentration 0(x, i), energy density and 
nutrient fields. 

To determine the contributions of activity to the 
stress, we need the forces associated with the ac- 
tive particles. Since there are no external forces 
on the system, the simplest active particle, on long 
timescales, is a permanent force dipole (see Fig. 
P). A collection of such particles [j] where the ath 
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FIG. 1: Force dipoles for (a) a rowboat, (b) a "bac- 
terium" with two fiagella [l4| . and (c) a motor on a fila- 
ment 

particle, centred at R Q , has point forces of strength 
/ and directions ±n a situated at R a + bn a and 
R„ + b'ii a , leads to a force density 

F a ~ -(6 + &')/V-^n a n Q 5(r-R Q ) 

a 

+ ^M/VV:^nAn^(r-R fl ) 

a 

= V-<x a , 

which defines the active stress cr a . For bacteria 
swimming at speed v$ in a fluid of viscosity 770 1 
/ ~ rjobvo. Both polar (6 ^ b') and apolar (6 = b') 
particles disturb the fluid; the former (a "mover") 
induces a nonzero fluid velocity at its centre and 
hence moves, the latter (a "shaker"), by symmetry, 
cannot. Note that <x a is insensitive, at lowest or- 
der in gradients, to the asymmetry b — b' : movers 
and shakers have the same far-field fluid flow, and 
an isotropic collection of either or both should have 
similar rheology. Since the force dipole determines 
an axis for each active particle, the natural defini- 
tion </>(r)Q(r) = X) a (nan Q - |l)5(r - R a ) (where 
I is the unit tensor) of a local nematic order pa- 
rameter or alignment tensor Q associated with the 
activity lets us explore the rheological consequences 
of spatiotemporal correlations in the activity by a 
simple generalisation of nematodynamics |15j . We 
have thus established that the active contribution 
to the deviatoric (traceless symmetric) stress |l6l 

er a — (l/3)\Trcr a = WQ + W 2 Q 2 + .... (2) 

where the constants W, W2 ~ (b + b')f<p characterise 
the strength of the elementary force dipoles, and the 
sign of W has vital rheological consequences. The 
relation J2J i s a t the heart of the novel mechani- 
cal properties of active systems [4(. Even without 
equations of motion, J3J tells us why an active sus- 
pension with long-range nematic order is different 
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FIG. 2: Discs (a) and (c) and rods (b) and (d) with ac- 
tive force densities attached along their symmetry axes, 
under shear (horizontal arrows). The parameter W > 
in (a) and (b) and < in (c) and (d). 

from its passive counterpart. Both have Q ^ 0; the 
passive nematic, bound by Pascal's Law since it is 
an equilibrium liquid despite its orientational order, 
has a purely isotropic mean stress, i.e., a pressure, 
whereas the active nematic has a nonzero mean de- 
viatoric stress, a truly nonequilibrium effect. 

Fyg. 121 shows what the parameter W means. In an 
imposed flow, in the absence of activity, discs (rods) 
tend to spend most of their time with symmetry axis 
along the compression (extension) axis of the flow 
[ri| . When activity is switched on, the flow induced 
by the intrinsic force dipoles will clearly oppose the 
imposed flow in cases (a) and (b), and enhance it in 
(c) and (d). 

For passive nematogens, Q is governed by a free- 
energy functional F[Q] containing polynomials in Q 
as well as Frank elastic terms ~ VQVQ, giving rise 
to a passive order-parameter stress 

cr OP = 3G G Q Q G (3) 

where G = -6F/SQ+(l/3)lTrSF/6Q is the nematic 
molecular field. The mean deviatoric passive stress 
J2J) is zero in both isotropic and nematic phases. 
For small nematic perturbations <5Q in the isotropic 
phase, F oc a J 4>Tr(8Q) 2 so that the stress fluctua- 
tion ~ a<j)5Q with a coefficient a which decreases on 
approaching the transition to the ordered phase. In 
active systems, the relation © between stress and 
order parameter does not arise from a free-energy 
functional, and the proportionality constant W has 
no reason to decrease with increasing nematic cor- 
relations. This difference will be seen to be crucial 
when we compare the pretransitional viscoelasticity 
of passive and active nematogenic suspensions. 

Including the viscous stress a" = —770 A + O(QVu), 



expressed in terms of the rate of deformation A = 
(l/2)[Vu + (Vu) T ] and the hydrodynamic velocity 
field u = g/p for a system of density p, the total 
deviatoric stress er in the active case can be written 



—a i —v 
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plus a noise source uncon- 



strained by a fluctuation-dissipation theorem since 
this is a nonequilibrium system. This defines com- 
pletely the equation of motion for the momentum 
density g. 



The coarse-grained equation of motion for u follows 
from that for Q which when linearised involves only 
ter ms |4| of a form present in passive nematodynam- 
ics Il5|7 
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where r is the activity correlation time, D is a diffu- 
sivity which in passive systems would be the ratio of 
a Frank constant to a viscosity, Ao is a "reversible" 
kinetic coefficient ^t|j f is a traceless, symmetric, 
spatiotemporally white tensor noise with variance 
Nq, representing thermal or active fluctuations, and 
the ellipsis includes the coupling of orientation to 
flow. 

We are now ready to calculate the linear viscoelastic 
properties of our active suspension. In the isotropic 
phase, Eqs. |J2J, © and (@J, linearised and ap- 
plied to spatially uniform oscillatory shear flow at 
frequency to in the xy plane, imply 
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which defines the storage and loss moduli G'(u) and 
G"(uj). This is the claimed active enhancement or 
reduction n act oc Wt of the effective viscosity at zero 
shear-rate and zero frequency. Activity enhances 
viscosity in Fig. El (a) and (b), since W > 0, and re- 
duces it in (c) and (d) (W < 0). For W > © tells 
us that the viscosity grows substantially as the sys- 
tem approaches a transition to orientational order 
(which is in general continuous for active vectorial 
order Q), i.e., as r is increased. By contrast, in a 
passive system approaching a nematic phase the ex- 
cess viscosity ~ ar is roughly constant since r oc 1/a. 

Eq. JSJ also predicts strong viscoelasticity as r in- 
creases. For passive systems W = 0. Since a oc r _1 , 
G'(u>t 3> 1) decreases as Xq7]o/t. There is little 
viscoelasticity near an equilibrium isotropic-nematic 
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transition. For active systems, by contrast, W is in- 
dependent of r and of proximity to the transition. 
Thus, as r grows, 



G'(ujt > 1) ~ W 



(6) 



independent of r and, of course, the dynamic range 
over which elastic behaviour is seen increases. At 
equilibrium, one would expect such strong viscoelas- 
tic behaviour from a fluid or suspension near trans- 
lational freezing, not near orientational ordering. 

The contribution W^Q to the deviatoric stress in ac- 
tive systems modifies sharply the stress vs rate flow 
curve. To see this, start with a passive sheared ne- 
matogenic system in the isotropic phase near 
the transition to a nematic. Qualitatively, as the 
shear-rate is increased from zero, Q increases ini- 
tially linearly, then more rapidly and then essentially 
linearly again, leading to shear-thinning [l9j |. If we 
switch on activity, with a positive value of W, the 
rapid increase in Q implies an equally rapid increase 
in <r a . This will at the very least mitigate the shear- 
thinning and, if strong enough, will lead to shear- 
thickening. Alternatively, a system with W < will 
enhance the unstable shear thinning. Note that the 
sign of W can be got from an independent experi- 
ment at low concentration, simply by seeing whether 
switching on activity increases or decreases the vis- 
cosity. Thus, the effect of activity on the zero fre- 
quency shear viscosity predicts the shear-thickening 
or -thinning nature of the active suspension. 

We now calculate active fluctuation corrections 
to the shear viscosity. Eq. J5J contributes 
an active force density ~ W2VQQ to the mo- 
mentum equation, whose effect on viscosities, 
at one-loop order, is of the form Ari(u) ~ 

W^/d 3 fcd£exp(^)G Q (M)C Q (-fc,i), where G Q 
and Cq are respectively the propagator and corre- 
lation function of Q. From 
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for wr> L 



(7) 



Expressing the noise strength in terms of an ef- 
fective temperature T e ff, and assuming on dimen- 
sional grounds r/ /W ~ t , Nq ~ kBT e ff/T , and 
D ~ ^ 2 /t , t = Vo£ 3 /kBT e ff is the rotational relax- 
ation time of a single active particle, and I its typi- 
cal size, we see from Q that Ari(u})/rj ~ (ljtq)^ 1 ^ 2 , 
i.e., G*{uj) — \fiuj. 

All of the above effects are likely to be greatly en- 
hanced if the transition is to a polar-ordered phase, 




FIG. 3: Flow fields due to polymerisation and depoly- 
merisation at two ends of an aggregate 



since such a transition is expected .2] to be continu- 
ous, so that t can increase without bound. Further- 
more, since the bare timescale To is of order seconds 
for bacteria, the effects can be observed over a large 
dynamic range. 

Activity greatly enhances the noise temperature: on 
dimensional grounds the variance of cr a (k — 0, w — 
0) is ~ W 2 £ 3 t, with W ~ V u o/£, f° r active parti- 
cles moving with typical speed uq, correlated over 
a scale £ and time r. Equating this to A^T^ 77 
and estimating 77 ~ rj water — 0.01 poise, uq to be 
a bacterial swimming speed ~ 20/im/s, r ~ 1 sec 
(an E. coli run time) gives us an noise temperature 
T eff ~ 10 5 - 10 6 K, consistent with 20]. This will 
mean a thousandfold enhancement of the t~ d / 2 long- 
time tails in the autocorrelation of tagged-particle 
velocities. On timescales shorter than r, effects as- 
sociated with spatiotemporal correlations in Q [2lJ 
intervene. For a drop (23 or a film [^J of size L the 
tails will be cut off on a scale r v ~ pL 2 lit 2 r\. In |2(il |. 

T ~ T„. 

Rheology enters biology crucially through the active 
order parameter Q in several motility experiments 
which we discuss below. In gels imposed strains as 
well as elastic anisotropies enter in exactly the 
same way as A does in a fluid medium. This provides 
a natural explanation for elastotaxis, the ability of 
individual motile rod-shaped bacteria such as Myx- 
ococcus xanthus to orient with their long axes along 
the extension axis of an imposed elastic stress in 
their substrate [ill ], as well as active mechanosens- 
ing ^3] j where cells orient along the axis of greatest 
rigidity of an ambient gel. Cytoplasmic streaming 
|lCj . associated with the crawling of amoebae, arises 
naturally in our model. Fig. shows that an aggre- 
gate actively polymerising at one end and depoly- 
merising at the other has induced flow fields with ex- 
tensional and compressional axes interchanged. The 
resulting gradient in the active stress can be seen, 
in Fig. |31 to generate a mass flux from left to right. 
The effect will be enhanced by the fact that the de- 



polymerising end, with negative W, is shear thinning 
and hence more fluid. These arguments suggest why 
such streaming always accompanies amoeboid loco- 
motion. 

To summarise, we have constructed the general 
equations governing the rheology of suspensions of 
active particles, and derived several novel predic- 
tions, quantitative and qualitative. Our description 
is universal: only the values of parameters such as 
W, t and Ao distinguish the rheologies of a bac- 
terial suspension and a motor- microtubule extract. 
We look forward to tests of these predictions in ex- 
periments on living, reconstituted, or artificial [2^ ] 
active-particle systems. 
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